Quantum spin pumping at fractionally quantized magnetization state for a system 

with competing exchange interactions. 
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We study the quantum spin pumping of an antiferromagnetic spin- 1/2 chain with competing 
exchange interactions. We show that spatially periodic potential modulated in space and time acts 
as a quantum spin pump. In our model system, an applied electric field causes a spin gap to its 
critical ground state by introducing bond-alternation exchange interactions. We study quantum 
spin pumping at different quantized magnetization states and also explain physically the presence 
and absence of quantum spin pumping at different fractionally quantized magnetization states. 



1. INTRODUCTION 

An adiabatic quantum pump is a device that generates 
a dc current by a cyclic variation of some system param- 
eters, the variation being slow enough so that the system 
remains close to the ground state throughout the pump- 
ing cycle. The pumping physics gets more attraction 
after the pioneering work of Thouless [2, Quantum 
adiabatic pumping physics is not only related to the spin 
system but also related to the other systems like open 
quantum dots 0, 0, HI , superconducting quantum wires 
[1, 0], the Luttinger quantum wire [§] and also to the 
interacting quantum wire 

The motivation of our study is the following: we have 
understood from the previous paragraph that adiabatic 
quantum pumping may arise in different systems due 
to the presence of different pumping sources. Here we 
would like to study the adiabatic quantum pumping of 
a system that has not covered in any one of the pre- 
vious studies. We consider an antiferromagnetic spin- 
1/2 chain with competing exchange interactions to study 
the adiabatic quantum spin pumping. We consider both 
nearest-neighbor (NN) and next-nearest-neighbor (NNN) 
exchange interaction in the spin chain and only consider 
the presence of electric field that induces time dependent 
dimerization in both the exchange interaction. We also 
study the spin pumping at different states of magneti- 
zation. Our approach is completely analytical. We use 
Abelian bosonization and one-loop RG calculation to ex- 
plore spin pumpingphysics of this model Hamiltonians 
system. Shindou [lfj has studied only Heisenberg XXZ 
spin chain with NN dimerization. There are no compet- 
ing exchange interactions and also the effect of differ- 
ent states of magnetization on adiabatic spin pumping 
physics is absent. Shindou fl(]| has considered two per- 
turbations which opens a gap in the excitation spectrum. 
One of them is the bond-alternation exchange interaction 
which leads to the dimcrizcd state and the other one is 
the staggered magnetic field which locks the spin into a 
Neel ordered state. In his model, applied cyclic electric 
and magnetic field control staggered component of ex- 



change interaction and staggered magnetic field respec- 
tively. In our case there is only one perturbation which 
opens gap in the excitation spectrum. A par t of our 



model has some experimental relevance llj, Il2| . Sup- 
pose we have a spin-1/2 chain (like Cu-benzoate and the 
charge order phase of Y&4AS3) with unit cell containing 
two crystallographic inequivalent sites, where both the 
translation symmetry ( T Sj — » Sj+i ) and the bond 
centered inversion symmetry ( Ibond S%-j <-> Si+j+i ) are 
crystallographically broken. But the NN exchange inter- 
action J in these spin chain does not have any alternating 
component because the system has site-centered inver- 
sion symmetry which exchanges the NN bonds. If we 
apply electric field in a particular direction of the system 
[1 31 ] . then we may break the site-centered inversion sym- 
metry and it yields the bond alternation component in 
the NN exchange interaction. The additional interactions 
like, NNN exchange interaction and it's alternating com- 
ponents in the Hamiltonians are completely theoretical. 
We consider these terms in the Hamiltonians to study 
the nontrivial and interesting effects of these terms over 
the basic interactions. In sec II. we present the model 
Hamiltonians and general derivations. Different subsec- 
tions are for the different states of magnetization. Sec. Ill 
is devoted for conclusions. 



2. MODEL HAMILTONIANS AND 
CONTINUUM FIELD THEORETICAL STUDY: 



In our model Hamiltonians, we consider the presence of 
time dependent bond-alternation (dimerization) in both 
NN and NNN exchange interactions. We assume that 
the time dependence of dimerization is restored by the 
applied alternating electric field. In this section we do the 
all calculations, different subsections are for the special 
limit of this general derivations. Our model Hamiltonians 
are the following. 



Ha — J\ 



(S n (x)S n+1 (x) + S n {y)S n+1 (y) , 
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+AS n (z)S n+1 (z)) + J 2 2^(1 - S(t)(-l) n )S n ■ S n+2 
n 

- gf x B Hj2Sn(z) (1) 



h b = Ji2(i-«y(t)(-i) 

n 

+A 5„(^)5 n+ i(z) 



(S n (x)S n+ i(x) + S n (y)Sr, 

+ J 2 



where n is the site index, x, y, and z are compo- 
nents of spin. J\ and Ji are the nearest-neighbor and 
next-nearest-neighbor exchange coupling between spins, 
Ji, Ji > 0, A is z component anisotropy of NN exchange 
interaction, S(t) is dimerization strength, which appears 
as a time dependent parameter in our Hamiltonians, H 
is the externally applied static magnetic field in the z 
direction. The staggered component of exchange inter- 
action is arising due to the broken site centered inversion 
symmetry under a electric field in a particular direction 
[10|. A site-centered inversation operation with the sign 
of elecrtic field reversed that requires, S must be an odd 
function of electric field [Io| . 

One can express spin chain systems to a spinless fermion 
systems through the application of Jordan- Wigner trans- 
formation. In Jordan- Wigner transformation the rela- 
tion between the spin and the electron creation and 
annihilation operators are 5* = ip^tpn — 1/2 , S~ = 
i> n eycp[mY, n jZ-oo n i\ > S n = i>n ex P[- i7r E"=-oo n i\ > 
where rij = ip]ipj is the fermion number at site j. 



Hai = - 



Ji 
2 



E 



l/2)(^ +1 ^ n+1 
-1/2). 



H 



A2 



E 



(^ +2 ^n + h.C.)(V'Li^ 



n+1 



E 



- 1/2) (^ +2 ^ 



n+2 



1/2), 
(3) 

-1/2) 
- 1/2X4) 



of Jordan- Wigner transformation for NNN exchange in- 
teractions. 

Similarly one can also recast the spin-chain systems 
with Ji dimerization into the spinless fermions. The 
Hamiltonians are converted as follows: Hbi = Hai, 
Hb2 = Ha2 and 

HB3 = -± 5(t)^2(-l) n (V4 +1 Vn+^n+l) 



+J 1 AS(t)J2(- 



1) 



n+1 



l/2)(^n-l/8) 



(2) . . n . 

Here our Hamiltonians are different from previously stud- 
ied dimerization problem. In Ref. [l4j and Ref. 15 1 
have studied intrinsic dimerization for frustrated spin- 
1/2 antiferromagnetic chain. In these studies, there is 
no explicit dimerization. Totsuka [l6| and Tonegawa 
have studied the Hb Hamiltonian only. These is no 
spin pumping physics in any one of the previous studies 
0, H EE UM- There are few other studies [H M, [H 
based on model Ha , but there is no spin-pumping physics 
in any one of these studies. So the current work is more 
wide and advance. Our approach is completely analyti- 
cal, i.e., we explain the basic understanding of spin pump- 
ing physics of our model system. Before we proceed fur- 
ther for continuum field theoretical study of these model 
Hamiltonians, we would like to explain the basic aspects 
of quantum spin pump of our model Hamiltonians: An 
adiabatic sliding motion of one dimensional potential, in 
gapped fermi surface (insulating state) , pumps an integer 
numbers of fermions per cycle. In our case the transport 
of Jordan- Wigner fermions (spinless fermions) is nothing 
but the transport of spin from one end of the chain to 
the other end because the number operator of spinless 
fermions is related with the z-component of spin density 



21|. We shall see that non-zero S(t) introduces the gap 



H A3 = - 



n 

j 2 s(t)j2(-i 



fh.c.XV^+iVVi 

1/2)^+2^+2 



at around the Fermi point and the system is in the in- 
sulating state (Peierls insulator). In this phase spinless 
fermions form the bonding orbital between the neighbor- 
ing sites, which yields a valance band in the momen- 
tum space. It is well known that the physical behav- 
ior of the system is identical at these two Fermi points. 
From the seminal paper of Berry [22T ]. One can anal- 
yse this double degeneracy point. It appears as source 
and sink vector fields defined in the generalized crys- 
tal momentum space (22[. B n {K) — V k X A n {K), and 
A n (K) = £ < n(K)\V K \n{K) >, where K = {k,5(t)). 
Here B n and A n are the fictitious magnetic field (flux) 
and vector potential of the nth Bloch band respectively. 
1 — l/2The degenerate points behave as a magnetic monopole in 
the generalized momentum space, whose magnetic unit 
— l/2(^) an ^ e snown t° be 1 [lo| . analytically 



There is a difference between the first term of Eq.3 with 
the first term of Eq.4 and Eq.5. This difference arises due 
to the presence of an extra factor e~ mn '+ 1 in the string 



dS-B± = ±1 



(7) 



Positive and negative signs of the above equation are re- 
spectively for the conduction and valance band. Conduc- 
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tion and valance bands meet at the degeneracy points. 
Si represent an arbitrary closed surface which enclose 
the degeneracy point. In the adiabatic process the pa- 
rameter S(t) is changed along a loop (T) enclosing the 
origin (minima of the system). It is well known in the 
literature of adiabatic quantum pumping physics that 
two independent parameters are needed to achieve the 
adiabatic quantum pumping in a system 0|. Here one 
may consider these two parameters as the real and imag- 
inary part of the fourier transform of dimerized potential. 
When the shape of the dimerized potential will change in 
time, then it amounts to change the phase and amplitude 
in time. The role of adiabatic parameters are not explicit 
in our study. We define the expression for spin current 
(J) from the analysis of Berry phase. Then according to 
the original idea of quantum adiabatic particle transport 

00, nam, the total number of spinless fermions (J) 

which are transported from one side of this system to the 
other is equal to the total flux of the valance band, which 
penetrates the 2D closed sphere (52) spanned by the T 
and Brillioun zone [icj |. 



dS ■ B +1 =1 



(8) 



We have already understood that quantized spinless 
fermion transport is equivalent to the spin transport 21 1 . 
We will interpret this equation more physically at the end 
(Eq. 18) of the next section. This quantization is topo- 
logically protected against the other perturbation as long 
as the gap along the loop remains finite 



1( 



In the following paragraph, we do the continuum field 
theoretical studies of spin pumping for different magne- 
tization states and explain the stabilization of quantized 
spin pumping against z-component of exchange inter- 
actions, and also from the intrinsic dimerization (when 
J2 > 0.241Ji [24j). We recast the spinless fermions op- 
erators in terms of field operators by this relation 



tp(x) = 



(9) 



where ipn{ x ) an d iPl(x) describe the second- quantized 
fields of right- and left-moving fermions respectively. In 
absence of magnetic field (H — 0), kp = ± tt/2, how- 
ever we are interested to study the systems in presence of 
static magnetic field. Therefore we keep Fermi momen- 
tum as arbitrary k F . One can simply absorb the finite 
magnetization in a shift of field <j) by 4> — <j) — nmx, 
where m =< S z > . In presence of magnetic field Fermi 
momentum and magnetization (m) are related by this re- 
lation, kp — ^(1 — 2m) [25J. We want to express the 
fcrmionic fields in terms of bosonic field by this relation 



1p r (x) 



-i (r<p(x) - 9(x)) 



(10) 



r is denoting the chirality of the fermionic fields, right (1) 
or left movers (-1). The operators U r are operators that 



commute with the bosonic field. U r of different species 
commute and U r of the same species anticommute. </> 
field corresponds to the quantum fluctuations (bosonic) 
of spin and 9 is the dual field of (p. They are related by 
this relation 4>r = 9 — tfi and 4>l = 9 + if). 

Using the standard machinery of continuum field the- 



ory 



251 ] . we finally get the bosonized Hamiltonians as 



Ho = v 



[ L *Jl{^ . H 2 
Jo 27T 1 



[d x ^x)f 



91 



dx : [d x <j) L {x)] 2 : + : [d x M x )Y 



+ || J dx [d x (j> L {x)){d x <f> R {x)) (11) 

Ho is the gapless Tomonoga-Luttinger liquid part of 
the Hamiltonian with vq — sinkp. The analyti- 
cal expressions for g\ and gi (related with the for- 
ward scattering of fermionic field) are the following. 
gi = 2(A — 2J2) sin 2 kp + 2 J2 sin2kp(n + sin2kp) 
g 2 = 4(A - 2J 2 ) sin 2 k F + 4 J 2 sin 2 2k F . 

Analytical expressions for different exchange interac- 
tions of Hamiltonian, Ha, are the following. 



Hj2Ci = „ V 9 / dx '■ cos[4:VK(j)<x) - (G-^k F )x - Ak F a] : 
2TT z a A ' 



(12) 



JiA 



H 



J2C3 



27r 2 a 2 



2ir 2 a 2 



Hj2C2 = 2" / dx '■ cos[4VK(j>(x) + {G—ik F )x - 2k F a] : 



(13) 



dx : cos[(-k— 4k F )x + 4v K <p(x) — 4k F a] : 

(14) 



Where G is the reciprocal lattice vector. Eq.12 and Eq.13 
are presenting the umklapp scattering term from the NN 
and NNN antiferromagnetic exchange interaction, Eq.14 
is appearing due to the presence of dimerized interaction. 
Similarly one can also find the expressions for Hb Hamil- 
tonian. Analytical expressions for K is the following. 



K 



1 — (8/tt) J2 sin 2 k F + AJ2 coskp 



1 + (4/ir)AsinkF + 4J 2 cos/c_f(1 + 2/irsin2kF) 

(15) 

vq and K are the two Luttinger liquid parameters. Dur- 
ing this derivation we have used the following relations: 
Pr/l = =^d x (j) R /L{x) and [<j>(x) , H(x)) =iS(x-x'), 
where n(a;) = —'V6(x), is the canonically conjugate 
momentum. We have also used the following equations, 
S z (x) = a [ p(x) + (-iy M(x) } . The bosonized ex- 
pressions for p and M are given by p(x) = — -j= d x <f>(x) , 

M(x) = cos(2(f>(x)) . Similarly one can calculate the 
analytical expressions for J\ dimerization. Here we have 
expressed our all expressions in terms of bare phase field 
(</>), by using the conventional practice of continuum 



1/2 



field theory [25J. During these derivations we assume 
that Ji 3> J2, S. J2 is in the unit of J\. Here we neglect 
the higher order of a than a 2 . 



2.1 CALCULATIONS AND RESULTS FOR M=0 
MAGNETIZATION STATES: 

At first we discuss m = magnetization state, it cor- 
responds kp = ± ir/2. Here we study both the effect 
of XXZ anisotropy (A) and the spin-Peierls dimerization 
{5(i)). The effective Hamiltonian for J2 dimerization be- 
come, 



H A = H + ( 



2tt 2 c 



) I dx : cos[AvK <j>(x)] 



(16) 

In this effective Hamiltonian (Eq.16), there is no con- 
tribution from dimerized interaction (kp = n/2 limit of 
Eq. 14), due to the oscillatory nature of the integrand 
(it leads to the vanishing contribution.). But the con- 
tribution of dimerized potential is present in the NN ex- 
change interaction. Similarly the effective Hamiltonian 
for Ji dimerization become, 



H B = H + ( 



JiA, 



27r 2 a 2 



) l dx: cos[AVK 4>{x)] 



+ / dx ■ cos[{2VK <j>{x)} : . (17) 

J 

This dimerization contribution for NN exchange interac- 
tion has originated from the XY interaction. This dimer- 
ization is the spontaneous dimerization, i.e., infinitesimal 
amount of 5(t) is sufficient to produce a gap around the 
Fermi points. The other two contribution of JI dimeriza- 
tion are from XXZ anisotropy of NN exchange interaction 
and z-component of NNN interaction. Fig. 1 shows the 
variation of K with ^ for different values of A. We 
observe from the figure that K is the function of A for 
fixed -j^ (inset of Fig. 1 shows the appearence of intrin- 
sic dimerization as a function of ^ and A. There are a 



few studies [26|, [27J on the phase seperation between the 
spin-fluid and dimer order phase of frustrated spin chain, 
but our Hamiltonians are different from them). The sec- 
ond term of Eq. 17 is irrelevant when K is greater than 
1/2. So in this parameter space only the time depen- 
dent dimerizing field (third term of Eq. 17) is relevant 
and lock the phase operator at <\> = + -^=. Now the 
locking potential slides adiabatically (here the cyclic elec- 
tric field that produces the dimerization) . Speed of the 
sliding potential is low enough such that system stays 
in same valley, i.e., there is no scope to jump onto the 
other valley. The system will acquire 2-7r phase during 
one complete cycle of external electric field around the 
loop encircling the minima of critical ground state, pro- 
duced by the dimerizing field. This is the basic mecha- 



m= 1/2 




0.3 
.4 



FIG. 1: Luttinger liquid parameter (K) versus -j 2 - for different 
magnetization state. A constant dashed line at if = 0.5 is for 
eye guide line, which separating the spin fluid and dimer or- 
der instability state. D-O: dimer order instability state, S-F: 
spin fluid Luttinger liquid phase. A: Here we only focus at 
m = magnetization state. Uppermost (dotted) curve is for 
A = 0.3 and the lowermost (solid) curve is for A = 1. The 
intermediate curves are for A = 0.4, 0.5, 0.6.0.8, 0.9 respec- 
tively from upper to lower one. Inset shows the separation 
between the spin fluid and dimer order instability state by 
a critical line. Here we present the shift of the (^) with 
A. B: Here we present the curves for different magnetiza- 
tion states, m — 1/2 is independent of A and ^ in contrast 
with m = 1/4, 0. The curves for m = 1/4, plateaus are for 
A = 0.4, 0.5, 0.6 respectively from upper to lower one. 



nism of spin pumping of our system. This expection is 
easily verified when we notice the physical meaning of the 
phase operator (<j) (x)). Since the spatial derivative of the 
phase operator corresponds to the z-component of spin 
density this phase operator is nothing but the minus of 
the spatial polarization of the z-component of spin, i.e., 
P s * = — 5Zj = i jSj Z . Shindou has shown explicitly the 
equivalence between these two consideration [lfj . During 
the adiabatic process < <f>t > changes monotonically and 
acquires - 2tt phase. In this process P s z increases by 1 
per cycle. We define it analytically as 



SP/ - J dP s z =~^J dxd x < <P{x) >= 1 



(18) 



This physics always hold as far as the system is locked 
by the sliding potential and A < 1 [10]. This equation 
(Eq. 18) for spin transport is physically consistent with 
the Eq. 7 (based on Berry phase analysis) of spin cur- 
rent. The quantized spin transport of this scenario can 
be generalized up to the value of A for which K is greater 
than 1/2 . In this limit, z-component of exchange inter- 
action and also the intrinsic dimerization has no effect on 
the spin pumping physics of applied electric field induced 
dimerized interaction of Hb Hamiltonian. 
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2.2 Calculations and Results for m — i 
Magnetization States: 

Here we discuss the physics of quantum spin pumping 
of a finite magnetization state. We are considering the 
magnetization state at m = 4, it corresponds kp — ±?. 
The effective Hamiltonian for Ji dimerization become, 

Ha = H ° ~ / dX '■ cos ^ ^ : 

dx : cos[^Vk (f>(x) + (3}: . (19) 



+ 



A 
2ir 2 a 2 



Where A 



J! 2 A 2 



J 2 \ P = tan- 1 ^. Ap- 
parently it appears from the general derivation of section 
(II), that the second and third terms of Eq. 19, will be 
absent due to the oscillatory nature of the integrand but 
this is not the case when one consider the dimerized lat- 
tice. In dimerized lattice, reciprocal lattice vector G will 
change from 2ir to tt due to the change of the size of the 
unit cell. It become more clear, if one write these terms 
as J dx : cos[(G - ikp)x + 4\/K 4>(x)] :. 

Similarly one can write the effective Hamiltonian for 
J\ dimerization: 

H B = H + (4=5§r) / dx ■ sin[4VK <j>(x)\ : 



2ir 2 a 2 



A 



2n 2 a 2 



dx : coslAVK 4>{x) +/?]:. (20) 



The analytical structure of Eq. 19 and Eq. 20 are the 
same, i.e., the coefficient of the field 4>{x) is the same 
for all sine-Gordon coupling terms. The renormalization 
group equations for these type of interactions are 2^, 28| . 



dK 
dlnL 



dd(t) 
dlnL 



AK)8{t) 



(21) 



(22) 



It appears from these RG equations that to get a relevant 
perturbation, K should be less than 1/2. It reveals from 
Fig. IB that K is exceeding the relevant value in our 
region of interest to mature criteria for spin pumping. 
So the dimerization strength should exceed some criti- 
cal value {5 C ) to initiate the spin pumping phase. These 
two equations are the Kosterlitz-Thousless equation for 
the system in this limit. At the critical point, system 
undergoes Kosterlitz-Thouless transition [2a, |28| • Since 
the system flows to the strong coupling (dimer-order) as 
the dimerization strength exceeds some critical value (5 C ) 
initially, we have to guess the physics of this phase. We 
analyze the system in the limit S — > ±oo and K — ► 0. 
In this limit all sine-Gordon couplings are relevant but 
the value of <fi is pinned at the minima of cos(4\^K(f>) for 



NN dimerization and of sin(4y/K(j)) for NNN dimeriza- 
tion because the dimerization strength is larger than the 
other couplings of the system, Hence it produces a deeper 
minima for the system. This parameter dependent tran- 
sition, from massless phase to massive phase, at T = 
is the quantum phase transition. This quantum phase 
transition occurs at the every magnetization state. So 
we conclude that the appearance of quantum spin pump- 
ing is not spontaneous like m = 0, rather dependent on 
the strength of the parameter. 



2.3 Calculations and Results for m = | and Others 
Fractionally Quantized Magnetization States: 



Now we discuss the saturation magnetization at 
rn = i (kp = 0). Kp = implies that the band is 
empty and the dispersion is not linear, so the validity 
of the continuum field theory is questionable. Values 
of the two Luttinger liquid parameters, vq and K, are 
and 1 respectively. It also implies that none of the 
sine-Gordon coupling terms become relevant in this 
parameter space. So there is no spin pumping for these 
fractionally quantized magnetization states. 



Here we present the explanation for the absence of 
other fractionally quantized magnetization state (like 
1/3, 1/5, 1/7 etc): A careful examination of Eq. 12 to Eq. 
14 reveals that to get a non oscillatory contribution from 
Hamiltonian one has to be satisfied Akp = G condition 
but this condition is not fulfilled for these fractionally 
quantized magnetization state. There are no sine-Gordon 
coupling terms. Hence there is no spin pumping physics 
for these fractionally quantized states of magnetization. 



3. CONCLUSIONS: 

We have presented the physics of quantum spin pump 
for different magnetization state of an antiferomagnetic 
spin-1/2 chain with competing exchange interactions 
along with bond-alternation interactions. Our study is 
completely analytical. We have observed that for some 
magnetization state spin-pumping is spontaneous and 
for some other it is not and also explain the physical 
reasons for the presence and absence of spin pumping 
for those states. 
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